Theory of Quantum Hall Nematics 
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Transport measurements on two dimensional electron systems in moderate magnetic fields suggest 
the existence of a spontaneously orientationally-ordered, compressible liquid state. We develop and 
analyze a microscopic theory of such a "quantum Hall nematic" (QHN) phase, predict the existence 
of a novel, highly anisotropic density-director mode, find that the T — long-range orientational 
order is unstable to weak disorder, and compute the tunneling into such a strongly correlated state. 
This microscopic approach is supported and complemented by a hydrodynamic model of the QHN, 
which, in the dissipationless limit, reproduces the modes of the microscopic model. 

PACS: 73.40.Hm, 73.50.Jt, 73.20.Mf, 64.70.Md 



Recent transport measurements on high mobiUty two 
dimensional (2d) electron systems Q have exhibited a 
striking anisotropy for Landau level (LL) filling factors 
close to v = N + 1/2, with > 4. A natural expla- 
nation of this anisotropy is the development of a local 
Id charge density wave (CDW) "stripe" state which ex- 
hibits spontaneous orientational order, that can be fur- 
ther pinned by weak crystal- or in-plane magnetic fields 
1^. Hartree-Fock studies had predicted the existence of 
inhomogeneous states with corroboration coming 

from exact diagonalization |^ and DMRG studies of 
small systems. Such inhomogeneous states of matter are 
ubiquitous in systems where there is a competition be- 
tween a repulsive, long-range (Coulomb) interaction and 
an attractive, short-range (exchange) interaction. 

While structurally similar to conventional CDWs in 
metals, in semiconductors the CDW is expected to be 
only weakly pinned to the underlying crystal because 
of the large disparity between the electronic and ionic 
densities. Therefore, in semiconductors a good starting 
point for the CDW state is the "quantum Hall smectic" 
(QHS) |p|-p^, in which translational and orientational 
orders are broken spontaneously, with crystalline fields 
acting as small perturbations. 

Although the QHS is an interesting state in its own 
right, an orientationally-ordered quantum Hall nematic 
(QHN) [pyr^] liquid is sufficient to account for the ob- 
served transport anisotropies. Additional motivation for 
studying the QHN is the instability of a 2d smectic to 
arbitrarily weak thermal fluctuations , which restore 
translational order and asymptotically reduce it to a ne- 
matic. Also, at T = sufficiently strong quantum fluctu- 
ations can unbind dislocations, driving quantum melting 
of the QHS into the QHN |]l5| and providing a continuous 
transition from the QHS to the isotropic, compressible 

= N + 1/2) electronic Uquid 

A conceptually useful picture of the QHN is a QHS in 
the presence of a plasma of free dislocations of density 
^J^. In conventional liquid crystals, where dislocations 
obey simple diffusive dynamics, this picture allows for a 



reliable derivation of the nematic hydrodynamics from 
that of the smectic However, the T = quantum 

dynamics of dislocations in the QHS is complex and is yet 
to be understood, as it depends upon the dislocation's 
charge, statistics, and the nature of the QHS state [ p5[ ; 
these subtleties are not addressed in treatments which 
use formal duality arguments g_8| or a Lorentz-invariant 
extension of the static model 

In this Letter we present a theory of the compress- 
ible quantum Hall nematic, which is guided by the chiral 
edge dynamics [^,|ri|,^ of the local smectic layers. The 
result is a novel quantum rotor model with a "soft" ki- 
netic energy, stemming from the compressibility of the 
QHN and the underlying non-inertial LL dynamics. Fur- 
ther support for the validity of our model comes from 
a complementary analysis of the finite T hydrodynamics 
of a charged nematic in a magnetic field, which in the 
collisionless limit reproduces the microscopic dynamics. 

In marked contrast to conventional liquid crystals 
where the constituent molecules are highly anisotropic, 
in the QHN it is the interaction-driven, dislocation-free 
mesoscopic regions of size £,d that play the role of "ne- 
matogens," with each region characterized by a local unit 
director field n(r) (normal to the local smectic layers). At 
scales longer than (which at T = 1 K we estimate to 
be of order 2000 A 0) the QHN Hamiltonian for the 
orientational degrees of freedom has the standard form 



H 



N 



+K,{V ■ hf + if3[fi X (V X n)]2 - (h • n)2}, (1) 

where Ki and if 3 are the Frank elastic constants pO[ , 
h is a symmetry-breaking field (e.g., the crystalline 
anisotropy or an in-plane magnetic field), x is the com- 
pressibility of this metallic state (including Coulomb in- 
teractions replaces by + STre^/g), i5/5(r) is the 
local electronic density operator, and in a single LL we 
have justifiably ignored the kinetic energy. 

The quantum dynamics of the problem is encoded in 
the commutation relation for the operators 5p and n. 
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This can be obtained by quantizing the classical dynam- 
ics of the edges of the local smectic layers (taken parallel 
to the X-axis). A fluctuation in the electron density 5p 
leads to a force —dx{x~^^P) parallel to the layers, which 
in the high field limit is balanced by the Lorentz force 
—eBii; a uniform compression along the stripe leads to a 
layer displacement uy, consistent with the fact that the 
LL-projected particle coordinates X and Y do not com- 
mute (equivalently encoded in the dynamics of chiral edge 
bosons P, pl] , p^ ) . Consequently, a nonuniform compres- 
sion along the stripe rotates the stripe through an angle 
9 = dxU determined by eB9 = d1{x~^5fi). Quantizing 
this dynamics on scales longer than |£9| leads to the 
commutation relation 

[(5/5(r), n(r')] = i^'z x n(r')a'<5(r - r'), (2) 

where I — ^J\/eB is the magnetic length and in our 
units h — ks = I- The corresponding QHN action Sn = 
J dtcPrCN is specified by a novel 0(2) "soft" quantum 
rotor model, with a Lagrangian density 

+Ki{V ■ n)2 + ^3 [n X (V X n)]' - (h • n)^}, (3) 

where L = L^z is the effective angular momentum 
operator conjugate to n, with commutation relation 
[^^(r), n(r')] = i('^z x fi(r')^(r — r') and defined by 
Sp — d^Lz- In addition to the first "Berry's phase" term 
associated with the conserved 2-component of the an- 
gular momentum, a notable feature of the QHN is that 
the kinetic energy (the second term) vanishes at long 
wavelengths, a property which can be traced to the non- 
inertial 2d electron dynamics in a strong magnetic field. 

Having defined the model, we can now assess the ef- 
fect of fluctuations on the QHN. Parameterizing the di- 
rector as n = (sin cos 0), we compute the (imaginary) 
time-ordered orientational correlation function C(r, t) = 
{%.9{r,T)9{O,0)), whose Fourier transform C'(q, w„) at 
Matsubara frequency aj„ = 2'imT can be measured in de- 
polarized dynamic light scattering experiments. We find 

C(q,c^„)=T^V^^, (4) 

with an anisotropic dispersion eq for orientational fluc- 
tuations about a y-directed nematic state given by 

eq = ±eX-^'\l^K^ql + K:,ql + h\ (5) 

exhibiting a line of nodes along the nematic order. In con- 
trast to naive expectations based upon experience with 
conventional nematics, in the QHN, because of the strong 
magnetic field, this orientational mode remains gapless 
(cq ~ even in the presence of an ordering field h — hy. 



For h = and in a single elastic constant approximation 
{Ki = K3 = K) the spectrum is oc ql\q\ (oc 
with Coulomb interactions). 

Standard analysis shows that at T = the equal- 
time connected correlation function Cc(r,0) = C(0, 0) — 
C(r,0) saturates at Cc(oo,0) = 9^^^^ - cie^/iQy/lt^) 
(ci — 0(1) constant), and therefore predicts that a 2d 
QHN exhibits true long-range order at T = 0, with 
the nematic order parameter ip2 — (e^*^) reduced by 
quantum fiuctuations from its T = classical value of 
1 to ^^2(0, ft.) ~ e-2«j^< 1. At low T, such that 
T <t: Tq = P-Q'^^jKlx, C'c(r,0) exhibits a plateau 

for Crf < r < = [P-jT)"''^ (Klxfl"^, but (for 

/i = 0) asymptotically crosses over to classical loga- 
rithmic growth ~ ln(r/^T) characteristic of quasi-long- 
range order of 2d classical nematics. For h ^ Q the 
order is long-ranged even at finite T, but is reduced 
from the quantum-renormalized T = value of '!/'2(0, h) 
down to V'2(T, ft) = i/'2(0,ft)i?(r,ft), with i?(T, ft) = 

exp [-S {^yao) f'/i8^'aoVK^)h Siy) = jf dzz{z + 
1)^1/2 J^2^rf5»cos2 6' [coth(fz\/7TTcos2 6i) - 1], A = 
£,ho/S,d, and ^/lO = K^^'^/h a strong-coupling orienta- 
tional pinning length associated with ft. In the quan- 
tum regime, T <^ Th = h^t"^ /{x^''^K) < Tq (equiva- 
lently, < ^0 < ^t), we find S{y) w C2/y^''^ (c2 ~ 

C(3/2)/(3^3/2) ^ Q_^g)^ ig^^i^g ^(y^ ^) ^ e-(T/To f'^^ 

with To = i^/^K^I^h/{cl'^x^'^)- Such nonanalytic 
(in T) thermal suppression of nematic order contrasts 
strongly with the naive classical spin-wave prediction of 
/mear reduction with T In the opposite, semi-classical 
Tft < T < Tq regime (equivalently, < Ct < 610), 
i?(T,ft) = /(^T/a), with f{y) « y^-'\ 7^2 = 2T/7:K, 
and = Ct(C°/Ct)^^^'*^''^^ the semi-classical pinning 
length, strongly thermally renormalized upward from its 
strong-coupling value ^^o- Even in this semi-classical 
regime, quantum effects are observable with T entering 
through both 772(1') and the quantum UV cutoff ^t(T), 
only crossing over to the conventional classical result for 
T > Tq, corresponding to ^t(T) < ^d- 

Next, we turn to the effects of quenched disorder on 
the QHN, with the dominant effects at long length scales 
coming from the random anisotropy field (5h(r) arising 
from, e.g., crystal defects and sample roughness, both 
randomly pinning the nematic director. Since the dis- 
order is static, in a linear "Larkin" approximation (valid 
at scales smaller than disorder persistence length ^a) the 
quantum dynamics drops out and we can adopt many re- 
cent results on the effects of disorder on classical liquid 
crystals ||2^ . For ft = we find that for d < 4 the nematic 
state is unstable to infinitesimally weak quenched disor- 
der; in 2d the local orientational order only persists out 
to the Larkin length « {2Tr)^^'^K/A^^'^, where A is the 
variance of the random orientational field Shx{r)5hy{r). 
On longer scales the full nonlinearity of the random field 
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must be taken into account; a preliminary analysis sug- 
gests the existence of a quantum nematic glass p9[ akin 
to its classical analog in nematic liquid crystals ||2^ . 

Tunneling can provide important spectroscopic infor- 
mation about the low-energy excitations. In order to 
calculate the tunneling current we need to construct an 
electron creation operator satisfying [Sp{r), 'i>^{r')] = 
^'t(r')(5(r — r'); the operator 



\i/t(r,i) 



V2^ 



i J, [D(r-T')e(r' ,t)+Arg(r~r')£^alL^{r' ,t)] 



(6) 

accomplishes this, with _D(r — r') a solution of £^d^D{r — 
r') = S{r — r') and a vortex attachment factor ensuring 
the fermionic statistics. The tunneling TV is determined 
by the Laplace-transform of the local imaginary-time 
Green's function G(r) = (r^«'(r, 0)^'t(r, r)) oc 6"*^, 
with $(r) = ^ J (CF(l-e-'ol^l)/eq. For short-range 

1/2 

interactions, the leading behavior is <&(t) « (|t|/to) 
(to = ^TT^K'^/'^x^^'^M'^U), independent of the ^r/L ratio 
(^3 = rfy^K^/h). TunneHng into the QHN is there- 
fore suppressed relative to both the isotropic v = 1/2 
Halperin-Lee-Read (HLR) state Q and the QHS ||ll|,|l2|. 
Coulomb interactions modify the large h limit, giving 
^(t) ^ |t|3/5, but leave the small h regime unchanged. 

Although above model of the QHN is quite appeal- 
ing and mathematically self-contained, on general physi- 
cal grounds a theory of a compressible QHN, based on a 
collective nematic orientational degree of freedom 9{r, t) 
alone is necessarily incomplete. This shortcoming is clear 
from the defining properties of the QHN state, a com- 
pressible anisotropic liquid^ in which free dislocations and 
their bound states of vacancies and interstitials (bubbles) 
constitute additional low-energy quasi-particle degrees of 
freedom. This contrasts strongly with the QHS, where a 
T = description in terms of collective phonon (chiral 
edge boson) degrees of freedom alone is physically con- 
sistent P , although not the only possibility jl^ . 

Consequently, we must extend above QHN model to in- 
clude these additional quasi-particle degrees of freedom. 
Having been unable to do so from first principles, we con- 
jecture that the quasi-particle sector of the full model 
is described by the HLR Lagrangian, £HLR[V'ja], with 
ip{r) and a(r) the composite fermion (CF) and statisti- 
cal gauge fields coupled to the nematic orientation 
collective degree of freedom n(r) via 

Ant = + (.9s|V • n|2 +gf,|V X h\^)^P^ . 

(7) 

This coupling is dictated by symmetry, with the first 
term the anisotropic CF mass proportional to the ne- 
matic order parameter Qij — riiUj — (l/2)(5.y [with 'ip2 = 
2(Qx£c + iQxy)]i the second and third terms the coupling 
of the CF density to splay and bend, and D = V — la. 



The success of the HLR state at = 1/2 |16| suggests 
that it might also be a good starting point at higher 
partially filled LLs, ly « 4- 1/2. Spontaneous devel- 
opment of orientational order (facilitated by screening 
of the Coulomb interaction by lower LLs) then leads to 
an orientationally-ordered HLR state. Such a state natu- 
rally coincides with the above picture of the compressible 
QHN as a melted QHS. The relation dxPi-v — dyO — 6(i(r) 
between the interstitial- vacancy density, = pi — py 
and dislocation density bd(jc) is reminiscent of the relation 
V X a = pcf . This intriguingly suggests the identifi- 
cation of the gauge field with vacancies and interstitials 
[i.e., a(r) = {0, pi-y)], and CF fermions with dislocations, 
lending further support to the above conjecture. 

The presence of these HLR quasi-particles will un- 
doubtedly lead to Landau-like damping of the nematic 
director dynamics. Nematic order will in turn induce a 
spontaneous Fermi surface distortion [|2^ and the con- 
comitant highly anisotropic transport observed in exper- 
iments [0,11 . The associated director n fluctuations will 
also contribute to the CF self-energies, modifying them 
relative to those of the isotropic HLR state. The quasi- 
particle sector of the model will allow for a second tunnel- 
ing channel, which, to lowest order, should add in paral- 
lel to the tunneling results discussed above. We leave the 
analysis of these and other related problems to a future 
publication | p^ . 

We now turn to hydrodynamics, which provides further 
support for our microscopic model of the QHN. Focusing 
on the Goldstone director mode n and conserved particle 
(p) and momentum (pv) densities (ignoring for simplic- 
ity the energy density), standard methods lead to a set 
of hydrodynamic equations for an orientationally-ordered 
2d liquid of charged rods in the presence of a magnetic 
field B — Bz, which to linear order are 



dtp = -poV • V, 



(8) 



-a 



N 



V 6p 



LUcVy 



-IxxVx - IxyVy + viidlv^ + lyssdyVx + {vi2 + iy33)dxdyVy, 
^ 2^0^' + (^) (^) - "^"^ 

-lyxVx - lyyVy + '^22dyVy + V^^dlvy + {vi2 + l'3.3)dxdyV^, 

dtO = -^d^vy —dyVa; - 7e-^, 

where i/jv is given by Eq. (^, lOc = eB/m, Vij's are 
viscosities, and and jij = "fo^ij + iQQij are the 
orientational and translational (momentum) relaxation 
rates. The reversible coupling A (determined by the mi- 
croscopic physics and generically expected to 
approach 1~ near the QHN-QHS transition), encodes 
the competition between the director's rotational con- 
vection and shear alignment. As in conventional nemat- 
ics, for |A| > 1 the director will shear-align at an an- 
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gle to a prescribed uniform shear flow v = yx, with 
tan^A = \/(A + 1)/(A — 1); for |A| < 1, such a steady 
state is impossible and fi will tumble anisotropically (via 
instantons). 

Analysis of Eqs. (^) in the coUisionless and Galilean- 
invariant limit {vij = "/e = lij = 0) leads to four modes, 
only two of which (corresponding to a mixture of den- 
sity p and nematic orientation 9) are hydrodynamic and 
which propagate with a dispersion 



2^1/2 



[(1 + \)ql + (1 - \)ql] ^K,q^^ + Ksql + h 



(9) 



Analysis of the eigenmodes reveals that this novel g'^ 
mode corresponds to clockwise/counter-clockwise direc- 
tor oscillations beating against density fluctuations. For 
A = 1 this mode is identical to the microscopic, T = 
dispersion in Eq. (||). This, together with the observa- 
tion that the long time limit of hydrodynamic equations 
can be obtained from a Poisson bracket between p and 
6*, generahzed to an arbitrary A suggests that the 
microscopic commutation relation, Eq. (|^), must be sim- 
ilarly generalized to 

[5p(r), n(r')] = "-fz x n(r') [(1 + A)^,^ + (i _ A)^,^] 5{r - r'), 

(10) 

with a corresponding replacement in the Lagrangian (|^) 
of d^L \ [(1 + \)dl + (1 - A)a2] L. A notable feature 
of the dispersion (|^) in the alignable |A| > 1 regime are 
lines of nodes along qy/qx = itani^;^, a property that 
should reflect itself in a variety of experiments. For in- 
stance, our tunneling result for A = 1 and short range in- 
teractions, <i>(T) ~ will become '-^ for ft, = 
and A 7^ 1, and for h ^ will become ~ logr (|A| < 1) 
and ~ logV (|A| > 1). 

In the presence of dissipation, at sufficiently long scales 
the two density-director modes become overdamped. For 
the Galilean-invariant case (7.^ = 0) the relaxation rates 



are w_ 



'iD_q and 



^iD+iK,ql 



K,qt 



with D± model dependent diffusion coefficients. How- 
ever, because QHN exhibits finite conductivity, at suf- 
ficiently long scales momentum relaxation becomes im- 
portant and the density and director dynamics decou- 
ple, leading to diffusive relaxation rates ujp — iDpq^ and 
LUg = iDg{Kiql + Ksq^). 

Although the velocity modes are non-hydrodynamic in 
the presence of a magnetic field, their relaxation deter- 
mines the linear longitudinal transport. In the isotropic, 
-02 = (HLR) state, this leads to the Drude model with 
an isotropic resistivity. In the QHN, "02 7^ and we find a 
resistivity anisotropy {pj;^ - Pyy)/{pxx + Pyy) = ip2lQho 
proportional to the nematic order parameter , making 
the h and T dependence of the order parameter experi- 
mentally accessible. 



In summary, by combining detailed microscopic and 
hydrodynamic analyses, we have proposed a model for a 
QHN, and argued that this novel orientationally-ordered 
electronic liquid crystal state may be realized in nearly 
half-filled high LLs. We predict the existence of a novel 
director-density mode, with a highly anisotropic q^ dis- 
persion, which remains gapless even in the presence of 
an ordering field. The nematic order is long-ranged at 
T = 0, but is destroyed by thermal fluctuations and 
quenched disorder. Clearly, many interesting theoreti- 
cal and experimental questions remain and now can be 
addressed in detail using the model presented here . 
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